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1. Introduction 

Let A be a random variable defined on a probability space (H, A, P), and let F be its distribution function 
with endpoint 

xp = sup{x S K, F{x) < 1}, 

and let E = 1 — F its survival function. 


Throughout the paper we suppose that E|A| < oo. The mean excess function e{u) of X is defined by (see, 
e.g., Kotz and Shanbhag [6], Hall and Wellner [5], Guess and Proschan |4j) 

[ [ F{t)dt if F{u) > 0 

(1.1) e{u) = ^{X-u/X>u)= I ^A)Ju 

y 0 whenever F(u) = 0. 


A natural way to estimate the mean excess function e(u) is achieved by using the plug-in method, that is 
replacing the survival function in (1.1) by its empirical counterpart, as did Yang [3]. 


Now consider a sequence Ai, A 2 ,... of independent copies of A. The plug-in estimator of e(u), for n > 1, is 


( 1 . 2 ) 


M) = 




T,7=l ^^hXi>u] 

\Xi>u\ 


where I[js:>u] = 1 if X > w and 0 otherwise. 
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MEAN EXCESS FUNCTION 


For notation convenience, we denote 


IPx(/«) = J fuix)dF{x) = J xdF{x) 


and 


IPa:(5«) = J gu{x)dF{x) = J dF{x) = F{u) 
where fu{x) = xl\x>u], gu{x) = and Pjc is the probability law of X. 


We also denote by P„ the empirical measure associated with the sample Xi, - ■ ■ , We have 

^71 1 ^ 

Vnifu) =and P„(5 k) = - X! ^[^">“1- 


2=1 


2=1 


Formulae (O and o» lead to 


Vxifu) 


— u if u < xf 


e{u) = { Vx{gu) 

0 if M > a;;’. 


and 


where X^ n = max Xi. 

l<i<n 


V’uifu) 

6rt('a) = i Pn(<?M) 


— U if U < Xn. 
if U > Xn 


One of the most important motivation of the study of the mean excess function comes from extreme value 
theory (EVT). Indeed this function e{u) is linear in the threshold u when E is a Generalized Pareto distribution 
(QPD) and this is quite a powerful graphical test for such distributions. 

By using the Vapnik-Chervonenkis (VC) classes and the entropy numbers technics, we have been able to 
establish that the empirical mean excess function e„(u) converges almost surely and uniformly. We showed 
that for any ui less than the upper endpoint of the distribution F, 

sup |e„(M) — e(it)| —)■ 0 a.s as n —)• oo. 

i2<ni 

Next, by using the modern theory of functional empirical process mainly exposed in m, we proved that the 
empirical mean excess function en(u) also weakly converges, that is 

{V^ieniu) - e{u)), It G /} {(G(/i„), u G /} . 

where G is a Gaussian process and G /} is a function family to be both precised later. 


Furthermore, using Talagrand’s inequality (see 0), and Mason and al. technics (see 0), we arrived at 
finding our best achievement: that is finding consistency bands for the mean excess function e(u). Precisely 
we establish that for any interval I = [uoi^i], ui being less than the upper endpoint of F and for any e > 0, 
we have for n large 

p(e„(u) - ^ < e{u) < e„(u) + u G /) > 1 - e, 

\ y/n y/n / 

where (E„)(„>i) is a non-random sequence of real numbers precised in Theorem and where F satisfies a 
very slight condition. 


These results allowed us to set graphical goodness of fitting test based on the empirical mean excess function 
and to apply this test to Dow Jones data. We found that the Generalized hyperbolic family distribution 
reveals, itself, to generally fit financial data. 
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In this remainder of the text, we are going to detail this outlined results, to demonstrate them, to make 
simulations studies about them, and finaly to apply them to financial data. 


The paper is organized as follows. We state uniform almost sure (a.s) convergence results in Section]^ and 
finite-distribution and functional normality theorems in Section 

SECTiON|^is devoted to setting a.s consistency bands for the mean excess function. In Section[^ simulation 
studies and data driven applications using Dow Jones data are provided. We finish the paper by a concluding 
section. 


Before we go any further, it is worth mentioning that, in the sequel, all the suprema, taken over u < ui, are 
measurable since the functions of u that we consider below, are left or right continuous. This means that we 
are in the pointwise-measurability scheme. Thus, even when we use the results and concepts in m, we do 
not need exterior either interior integrals or convergence in outer probability. 

2. Almost Sure Convergence 

In this section we are going to prove the uniform almost sure convergence of the empirical mean excess 
function by using Vapnik-Chervonenkis (VC) classes and bracketing numbers. 

Theorem 1. Suppose that E|g„(Ai)| < oo and E|/„(Ai)| < oo, then 

sup |P„(g„) - ^x{gu)\ -t 0 a.s as n-)■ oo 

U<Xf 

and 

sup |P„(/u) - Pv(/u)|-t 0 a.s as n^oo. 

U<.Xf 

For any fixed Ui < xp, for EA^ < -|-oo 

sup |e„(it) — e(u)| —>■ 0 a.s as n —)■ oo. 

U<Ui 


Proof. We observe that Fi = {gu,u < a;;’} is a class of monotone real functions with values in [0,1]. By 
Theorem 2.7.5 in m, the bracketing number Lr{Q)) is finite (bounded by exp(A/e), for every 

probability measure Q, any real r > 1, and a constant K that only depends on r). Since E| 5 u(Ali)| < oo for 
u < xp, is functional Glivenko-Cantelli class in the sense of Theorem 2.4.1 in uni, meaning that 

(2.1) sup |P„( 5 „) - Px(g«)|-t 0 a.s as n ^ oo. 

U<.Xf 


The class A 2 = {fu,u G [uo,ui]}, with ui < xp, is a Vapnik-Chervonenkis class with index ¥{^ 2 ) = 3 and 
its envelop is G = max(|/ii(,(a;)|, |/ui(a;)|). Then it satisfies the uniform entropy condition 2.4.1 in [10]. Then 
F 2 is a Donsker class and hence it is a Glivenko Gantelli class, that is 


( 2 . 2 ) 


sup |P„(/u) - Pv(/u)| ^ 0 a.s as n^oo. 

U<Xf 


To finish, fix ui < xp. Then for u < ui and n large enough, we have 


e„(u) - e{u) 


Pn(/«) _ Px(/«) 

Pn(5«) Px(5«) 

Pn(/«) Px(/«) Px(/«) Px(/«) 
Pn(5K) Pn(ff«) Pn(ff«) Px(5k) 


{Vu{gu))-\Vn{fu) - Px(/„)) - Px(/.) X 


Pn(gu) - Px(g«) 
Pn(gK)Px(g«) 
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Then 

(2.3) 


\en{u) 


e{u)\ < |P„( 5 „)|-' X |P„(/„) -Px(/J| + |Px(/„)| X 


|Pn(g«) - Px(g«)| 

\Vn{gu)Vx{gu)\ 


Let 


(2.4) 


e„ = sup |Pri(/„) - Px(/«)| and 6„ 

uKXfp 


From (2.1) and (2.2) above, we have 


sup |P„(5u) - Px(5'«)|- 

uKXjr 


Cji —^ 0 ci.s s^iid Sfi — y 0 cL.s , cls tl —y oc. 


Now for u < Ui , we have Px(5m) > Px(ff«i) and from (2.4) 


-Sn < Pn(5«) - Px(5u) < Sn 

-Sn + Px(5u) < Pn(5u) < <5„ + Px{gu), 

since Pn( 5 u) > Px( 5 u) - <5™ > 0 for n large enough, then (P„(g„))“i < (Px(g«i) -Sn)~^. 
We also have 


(2.5) 

Thus 

and then 


"x(/.)| = 


/ xdF{x) < / xdF{x) < / |a;|(iF(ai) = E|X| = a < oo. 

J u R «/R 


sup |e„(u) - e(u)| < e„ Px(5ui) - 

U<Ui *- 


sup |e„(u) — e(u)| —?► 0 a.s as n —>■ oo. 

U<U\ 


Px(5ui)(Px(5ui) - 5n) 


□ 


3. Asymptotic normality of en { u ) 

In this section, we are concerned with weak laws of the empirical mean excess process as a stochastic process. 
Hereafter {G’{g),g € G} denotes a Gaussian centered functional stochastic process with variance-covariance 
function 

r(5i,ff2) = J {giix) -Egi{Xi)){g2ix) -Eg2iXi))dF{x). 

Theorem 2. Let Xi, X 2 , ■ ■ ■ he iid rv’s with common finite second moment. 

Put I — [uojMi], with uq < ui < xp and define the functions o/t G K, 

Kit) =Ex{gu)~^fu{t) -'^x{fu)^x^{gu)gu{t) for u&L 

Suppose that F is continuous and satisfies 

f PiK - Fiv - S) 
limsup sup I ^ — “ — ——- 

Then the functional empirical processes {<G„((;„),u € /} and {G„(/u),u G /} weakly converge respectively to 
{G( 5 „), u G /} and {<G(/„), u G /} in £“(/). 

And {^/nien{u) — e{u)),u G /} weakly converges to {G(/i„),u G /} . 



Before we give the proof, we need this lemma. 
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Lemma 1. Let g be a finite measurable function defined on M such that < oo . Let uq < ui < xp- 

Define for any fixed S K and i5 > 0 


Let for a fixed n> 1, uG 


If sup 

uqKvKui ^ 


a'^{v,6)=[ [g{x)-'¥.{g{x)f dF{x). 

J v — S 

1 "" 

Sn{u) = [9(^3)\Xj>u) - ■ 


0 as 5 —>■ 0 and sup \g{x) — Eg(X)| < cxd, 

uq <x<ni 


then 


lim sup sup -P{ sup \Sn{u) — S'„(u)| > g) = 0. 

uq<v<ui n'>l 5 v — S<u<v 


■0 


Proof of Lemma [T| We fix u G M and consider a = sup |>S'„(u) — 5'„(u)|. Observe that for u < v, 

v—5<u<.v 


Sn{u) - Sniv) = ^ ^ {g{Xj)\u,v\{Xj) - E5(Xj)I]„,„](Xj)} . 
Since for all {u,v) G we have 

1 " 

\S^iv) - 5„(u)| < ^ E + |E5(X,)|] < oo, 

it comes that a is finite. So for any ^ > 0, we can find u G [v — S^v[ such that, 


(3.1) 


|S'„(u) - S'„(u)| > a - e. 


Now, let (5 > 0. Define for any p > 1, and consider Uj{p) = Uj = v — S + jS/p, j = 0, ...,p. 
Let us prove that for £ > 0, 

lim max |S'„(ii,) — iS'„(u)| > a — e. 

p—)-oo 0<j<p 


For each p > 1, let j such that 


Uj-l{p) <u< Uj{p). 


We have. 


|5'„(uj) - S'„(u)| > |S'„(m) - S'„(u)| 


> \Sn{u) - Sn{v)\ 


1 X V 

/— 'y '^, {9{^i)^]u,Ui(p)]iXi) (Wi) 

Rjip), 


by denoting 


We get from (3.1) 


Rj{p) 



y ( i9{^i)^]u,Uj (p)] i^i) ^9i^i)^]u,Uj (p)] 


1 


max — S'„(u)| > a — e + RAp)- 

0<j<p 


For a fixed n > 1, Rj(p) —)■ 0 as p —)■ oo, since the sequence of intervals (]it, Uj(p)])p>i decreases to the empty 
set as p —>■ oo. 
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Next, consider the collection points {uj{£), 0 < j < p,l < £ < p} and denote the set of its distinct values 
between them as {uj, I < j < rn{p)} . We still have \uj — Uj-i\ < 5/p. And we surely have for any £ > 0 

lim max ISniui) — Sn(v)\ > a — e 
p-S-oo 0<j<m(p) 

and then 

lim max — Sn(v)\ > a 

p-)-oo 0<j<m(p) 

and finally 

sup max \Sn{uj) — <S'„(u)| = a. 

p>l 0<j<m{p) 

By construction, max |S'„(uj) — iS'„(u)| is non decreasing in p. So, by the Monotone Convergence Theorem, 

0 < j<m(p) 

for any fixed u > 0, for any p > 0, 

(3.2) P( sup |5'„(u) - S'n(u)| >p) = lim P( max \Sniuj) - 5'n(T)| > p). 

v-S<u<v l<i<m(p) 


Put «.](^.)) .h>l. 

i—1 

We have 


m(p) 

^ ^ Zfi — 
h=j 


with 


m(p) 

^/n{Sn{v - 5) - Sn{v)) = Y Zi = Tm{p) = T(n, u, 5). 

2=1 


We observe that {Ti, T 2 ,..., Tm(p )} are partial sums of i.i.d. centered random variables so that the form 
a submartingale. By the maximal inequality form submartingales, for any fixed p 


( max |S'„(uj) - S'„(u)| > ry) = P( inax \Tj\ > pyf£i) < ) 




Er(n, u, 5)^ 


Since the right hand does not depend on p, we get by (3.21 


^P( sup \Sn{u) - Sn{v)\ > p) < I .^ ET{n, u, 5)'^. 
S 1! —5<22<U ^ 

Notice that T(n,u,8) is a sum of n i.i.d centered random variables with variance 


pv 

ki{v,5) = a'^{v,5) = / {g{x) -E{g{x)f dF{x) 

*/ v — 6 


and fourth moment 


K2iv,5)=f {g{x)-E{g{x))'^ dF{x). 

Jv-S 

Simple computations give (see the appendix O for a simple proof of that) 


E (r(n, u, (5))"^ = nK 2 (v, 5) + 3n{n — 1)k^(p, S). 
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By putting these facts together, we arrive at 

ip( sup |S'„(u) - S'„(i;)| > 77 ) < r]~ 

^ v—6<iu<iv 


_4 / nK2(v, S) + 3n(n — <5) 

d-n? 

Vo n d 


1 - - 
n 


Remark that 


K2iv,6) 

sup —— < 

Uq<V<Ui ^ 

< 


sup |g(a;) — E 5 (X)|^ x (5 ^ x sup f dF{x) 
\uq<X<Ui / Uq<V<U\ Jv — 5 


sup |ci(a;) - E 5 (X)| x sup 


uq <rE<iii 


Uo^V‘^^ 

F{v) — F{v — (5) 


Uq<V<Ui 


We finally get 


lim sup sup -P( sup |5'„(u) — 5'„(u)| > ??) = 0 

^>4 § v — S<u<v 


whenever lim sup 

^^”70 uo<v<ui 


S 


= 0 and sup \g{x) — E( 7 (X)| < + 00 . 

Uq<X<Ui 


□ 


This achieves the proof of the lemma. 


Proof of Theorem [ 2 ] 


By Theorem 2.7.5 in m applied to and by the fact that F2 is a Vapnik-Chervonenkis class, condition 
(2.5.1) is satisfied for both Fi and F 2 thus Fi and F 2 are Donsker classes. 

This may be used in a simple manner to get 

(3.3) An = max(sup |<G„( 5 u)| ,sup |<G„(/u)|) = Op(l,/) as n -)■ 00 . 

uei uei 

Denote the functional empirical process for any real function g by 

1 " 

Gnig) = -Eg{X,)} . 

Remind that for any Donsker class Q, the functional stochastic process {'Gnig),g € G} converges in law to a 
Gaussian and centered stochastic process {G{g),g G G} whose variance-covariance function is 

r(gi,52) = J{gi{x) - Eg{Xi){g 2 ix) - Eg 2 iXi))dF{x). 


En{gu) = Exigu) 


Gn{gu) 

y/n 


Enifu) — Ex{fu) + 


Gnifu) 

y/n 


We have, as n —)■ 00 
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Thus 


Vn{en{u) - e{u)) = 


lifu) Pa(/«) 


i{ 9 u) Vx{gu) 

^Uu) Px(/«) , Px(/.) Px(/«) 


i{gu) Pn(g«) Pn(g«) Px(5k) 


1 


Pn(gu) 

1 

Pn(5u) L 
1 


_/ s V^(Pn(5u) - Px(5(w)) 

*^(Pn(/«) ~ Px(/«) j ~ Px(/u) I 

^ ^ Pr,,(5«)Px(5«) 




G, 


{fu 


We find 


(Pn(5«))' 


Px {gu) 


Pn(g«) L 
^nigu)'' 


Px {gu 

Vxifu) , 

Px(ff«) 




X 

D—1 / 


1 +P^^(g„) X n X G„(gu) 


-1 


= 1-Px^(ff«)^”- X G„(5 u)+P^^( g„) X 0(^71 X G„(5 u) 


Since is a Donsker class, then sup |G„(g„)| = ||G„||jfi = Op(l,/). So 

uG/ 


(Pn(5«)) ^ = Px^(5«) 1 - Px^(ff«) ^ Op(l,/) 

Let us remind that /i„ = Px(5«)“V« - ^x{fuWx^{gu)gu- Then, for m e /, we get 

v^(e„(M) - e(u)) = Gn(^fu - ^ Px^( 5 «) - Px^( 5 u) x x Op(l,/) 

= G„(/iu) +G„(/i„) X Px^ffu) X X Op(l,/). 


We finally have 

(3.4) ^/n{eniu) - e{u)) = Gn{hu)+Gnihu) X or{l,I). □ 


Lemma 2. 


The class J -3 = 



fu 

Px (gu) 


Px(/„) 

Pi(5.)^“ 


,M e /| 


is a Donsker Class. 


At this step, we want to prove that J^ 3 = {h^, uq < u < mi} is a Donsker Class. Since we obviously have, 
by the Central Limit Theorem, finite distribution convergence of {G,^(/i„),m G /} to the stochastic process 
{G(hu),u G /} in i°°{F^), we only need to prove the asymptotic tightness of {G„(/i„),u G /}. 


In view of Theorem in 1.5.7 in [10) . it is enough to prove that 

limsuplimsup -P( sup |Gn(/i„) — G„(/it,)| > 77 ) = 0. 

itG/ n—fCJO 5 v—S<u<v 

Here, we apply Lemma [T| for the nondecreasing mesurable function g{x) = x and g{x) = 1. 


In both cases, we inspect the assumptions of this lemma and see that if g{x) = x 
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we get g(x) < g{ui) = ui for any uq < x < ui and thus 

^ ua< 


sup 

Uq<V<U\ 


i i9{x)-^g{x)fdF{x)\ 

<v<ui 0 V Jy_s / 


< |ni —E(X)|"‘x sup ( 

Uq<V<Ui ^ 


F{v) — F{v — (5) 

71 


2 

—y 0 us S —y 0, 


and 


sup \g{x) - E 5 (X)| < |ui| + |E(X)| < oo. 

xGl 


If g{x) = 1, the result is obvious. 


We can apply Lemma and we will get, 


and 


lim suplimsup -P( sup |<G„(/„) - <G„(/„)| > rj) = 0 

uGl n—^oo S v — S<u<v 


lim suplimsup -P( sup — Gn{gv)\ > v) = 0- 

uGl n—^oo S v—S<u<v 


But by Theorem 8.3 of Billingsley [2], p.56, and by Theorem 2.2 in Lo [7], these two previous equalities entail, 
that 


and 


lim sup lim sup P( sup |Gn(/«) - <G„(/„)| > rj) = 0 

0 —)-0 n—)-oo \u-v\< 6 ,{u,v)gI‘^ 

lim suplimsupP( sup |<Gn(5u) — <Gn(<?«)| > v) = Q- 

'S-s-OuG/ n-).oo \u-v\<S,{u,v)eI^ 


Next, we use the following development for {u,v) G P 

Px(/.) 


hi, hi, — 


Px {gv 

Px(/«) Px(/.) 

Px(5«)^“ Px(5«)^'’ 


x{9u)[fu ^x{9v)[U 

^X {9n)fu — ^x {9v)fv—( 

^ ^ V 

a{u,v) — 


b{u, v) 


We get 


a{u, v) = P3^ {gu)fu - Px i9u)fv + Px {9u)fv - Px {9v)f, 

= ifu - fv) X P7(5«) + (p7(ff«) - ^7(5^;)) X /, 
fu-fv Px(5«) - Px(5'i,) 


^x{9u) Px(5k) X Px(ffi,) 


X fv 


Then 


Next 


|Gn(a(«,u))| < X |G„(/„ - /,)| + X G„(/„). 


Px(5«) 


Px(5«) X Px(ffi;) 


,, , Px(/.) ,, Px(/.) ,, , Px(/„) „ Px(/.) ,, 

Px(/«) , rPx(/.) Px(/.)- 


= i9u - gv) X 


Px(5“) ^Px(ff“) Px(ff'«)- 


X gv 


Next, 
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Px(A) Px(/.) 
Pi (5.) Pi (5.) 


Px(/„) Px(/„) , Px(/„) 


Pi(5«) Pi (5.) Pi (5.) 
(p^ ■ ^ 


Px(/.) 

Pi ( 5 .) 


+ (Px(A)-Px(/.))x^ 

/ (Px(gt;)+ Px(g«)) X (Px(g^) - Px(gu)) \ 
^ Pi(g«) X Pi(g„) / 

+ (Px(« - Px(W) X 


Also, 


IG.W,-))! < W^>‘IGx(9.-9.)I 


+ ^tox)+Pj;to.)|x|Px(9,:)-Px(!,.)) ,, 


Pi(gM) X P^(g„)| 


|Px(/«) - Px(/,;)| X |G„(g„)| X 


Pi(gG)' 


For {u,v) e let us use the bounds of Pj^^(gu), Pjf^(g„), and Px(/«)- 


We obtain < (F(mi)) i P^^( 5 „) < (i^(ui)) ^ 


and finally, from (2.5), we get |Px(/«)| < E|X|. 


Thus by using these bounds and ( |3.3[ ), it comes that 

sup \Gnihu - K)\ < Bi X sup |<Gn(/« “ /«)! + -82 X SUp |<Gn(gu “ 9v)\ 

\u—v\<5,{u,v)£l^ \u—v\<6,{u,v)gI^ \u—v\<5,{u,v)gI^ 

+ (-B 3 X sup |Px(g«) - Px(g«)| + xB 4 X sup |Px(/«) - Px(/«)| ) A„, 

\ \u—v\<5,(u,v)^I^ \u—v\<5,(u,v)gI^ J 

where 


'Si = (F(«i))-i; 

S 2 =E|X| X (F(wi))-^ 

< S3 = fK))-2((s(ui))-2 + 2FK)); 

S 4 = (F(ui))-i 


An = max(sup |<G„( 5 u)| ,sup |G„(/„)|). 

uei uei 


Now we observe that 

sup |Px(g^)-Px(g«)| = sup |F(u)-F(u)| 

\u—v\<6,{u,v)gI^ \u—v\<5,{u,v)gI^ 
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and 


sup \Vx{fu) -Vx{fv)\ < 

\u—v\<S,{u,v)gI'^ 


< 


sup 

\u—v\<S,{u^v)^I^ 


sup 

\u—v\<S,{u,v)^I^ 


poo poo 

/ tdF{t) — / tdF{t) 

J U j V 


tdF{t) 


< maxduol, luil) sup \F{u)-F{v)\. 

\u—v\<S,{u,v)£l^ 

These quantities go to zero whenever F is continuous and hence uniformly continuous in I. Putting all these 

lihu — ^d)| —^ 0 J —y 0. 


facts together and using (3.3) yield 


sup sup 

n>l \u—v\<S,{u,v)£l‘^ 


Finally F 3 is a Donsker class, thus sup \Gnihu)\ = 0]?{1,I) and we get from (|3.4|) that 

u^I 


>/n(e„(u) - e(u)) = + op(l, /). 


□ 


This completes the proof. 


Now we are going to concentrate on consistency bands for the mean excess function. 

4. Consistency bands 

Now, we may use the uniform bands of the functional empirical processes based on Talagrand’s inequality 
(see HI) and the new methods introduced by Mason and al. [9] to obtain consistency bands of the mean 
excess function as follows. 


Theorem 3. Let Xi, X 2 ,-- ■, he i.i.d random variables with finite second moment. Put I = [uo,ui], with 
— 00 < Uq < Ui < Xp- 

We suppose that F is continuous and satisfies 

(4.1) hmsup sup -p- 

s^o {v,v-S)ei^ \ VS 

Then for any e > 0, there exists uq such that for n > uq, 



with 

(4.2) En 

and where 


A and Ai are universal constants. 


E E \ 

-^ < e(u) < e„(u) H — u S / J > 1 — e, 

Vn V'^ ' 

_ 1 Hi^xElXK 

i^(Mi) - Hi/Vn V ^ F{ui) /’ 

( Di = 2 AAi \/log 2 + Al 
\ D 2 = AAi Ml \/log Ml + Al, 

[ Ml = max(2,max(|iio|, |ui|)) 


The proof of this theorem is rather technical so we postpone it in the appendix 
we also state the fundamental Talagrand’s inequality. 


subsection 


7.2.1 where 
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Remark: The validity condition (4.1) is quite very weak and is satisfied by most of the continuous usual 
distribution functions. Indeed if F is absolutely continuous with respect to the Lebesgue measure with 
derivative function /, we get by using the mean value theorem, 

2 

< 5 X sup f^ix). 

xG [v—5,1!] 

But sup f^(x) < oo whenever / is continuous, by a simple argument from real analysis. This allows 

xG [i! —5,?;] 

consistency bands for a huge number of absolutely continuous distribution functions. All the examples in the 
Section]^ are devoted to simulations satisfy (4.1) through this argument. 


F{v) — F{v — S) 

71 


Now, we are going to focus on the applications of our results. 


5. Simulations and applications 
5.1. Introduction. The Mean excess function can be used in two ways : 

• First, it can be used to distinguish heavy tailed models distribution and those with light tailed distribu¬ 
tion. An increasing mean excess function e{u) indicates a heavy-tailed distribution and a decreasing mean 
excess function e{u) indicates a light-tailed distribution. The exponential distribution has a constant mean 
excess function and is considered a medium-tailed distribution . 


Then the plot of the mean excess function tends to infinity for heavy-tailed distributions, decreases to zero 
for light-tailed distributions and remains constant for an exponential distribution. 


• Secondly, it can be used for tail estimation with the help of the generalized Pareto distribution which 
can model the tails of another distribution. 

Let Fu(x), the excess distribution over threshold it, defined by 

Fu{x) = P(A — u < x\X > It) 

with 0 < X < xp — u, where xi? < oo is the right endpoint of F. 


By using Theorem 7.20 in [I] , a natural approximation of F^ is a generalized Pareto distribution QPD{^,f3) 
which mean excess function is given by 


(5.1) 


3(u) = 


provided that ^ < 1. 


P , ^ 

1-c 1-7 

If the empirical mean excess function plot looks linear, we can fit a QPD(^, (5) model whose parameters can 
be estimated by means of linear least squares method : given data {(iti, yi),..., (itn, yn)}, where Ui = Xi and 
Vi = Gn{ui), 1 = 1,..., n, we estimate the parameters ^ and /3 to be 


where 


a -)- 1 


h 

a -|- 1 ’ 


with u = 


1 

n 


n 

Ui and y 

i^l 


n n n 

^ i—\ i—1 2=1 j r — ^— 

a = - - - - -^ and b = y- au, 

2 = 1 2 = 1 


1 


are the sample means of the observations on u and y, respectively. 


As far as we are concerned, our goal is to estimate the mean excess function by consistency bounds. 
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In the remainder of this section, we are backing on the empirical mean excess function {emef for short) to 
construct graphical tools goodness of fit test. 


In the first step we are considering a large set of distributions for which we draw the average emef. That 
means that we fix a distribution function and consider n = 6000 samples from it, each sample size is 4000. 
Next we compute the average of the n = 6000 empirical mean functions. 


The graphs of these average mean empirical functions would serve as stallions in the following sense: each 
other sample having an alike emef will suggest such an underlying distribution. 


We will use, as a special guest, the generalized hyperbolic {Gh for short) family of distribution functions. 
Nowadays, this family is very important in financial modeling. 


In a second step we will try to use the obtained graphs as stallions for real data. 

In this paper, we focus on monthly returns and log-returns of Dow Jones data. We will see that these data 
strongly suggest Gh model. 


This section, beyong financial data, shows how to use the emef for goodness of fit testing purposes. It opens 
a great verity of applications for differents types of data. 

5.2. Usual distributions. To assess the performance of our estimator, we present a simulation study. We 
draw simulated emefs for standard distributions and next for Gh family of distribution functions 

5.2.1. Emef for standard distributions. We consider some simple models that are listed in the TABLE[^below 
where the used parameters are specified and the emef figures corresponding to each choice are displayed. 


Distributions 

Parameters 

Figures 

GPD 

^ = 0.25, /3 = 1 

Figure 

f = -0.75, /3 = 1 

Pareto 

a = 7, A = 3 

Figure 

Exponential 

A = 2 

Weibull 

P = 1, T = 3.6 

Figure 

/3 = 1.5, T = 0.2 

Burr 

Gomberz 

a = 0.5, A = 0.5, T = 5 
a = I, A = 0.5 

Figure 

Gamma 

a = 2, /3 = 0.001 

Figure]^ 

Beta 

II 

II 

to 

Lognormal 

^ = 0, (7 = 1 

Figure]^ 

Normal 

^ = 0, (7 = 1 

Laplace 

^ = 0, CT = 1, r = 0.5 

Figure 

t Student 

u = 5, /r = 1 

Figure 

Gauchy 

II 

o' 

II 


Table 1. The emef for standard distributions 
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Figure 1. The emef for two generalized Pareto distributions : the left panel concerns the 
one with the parameters ^ = 0.25, /? = 1 and the right panel concerns the one with the 
parameters ^ = —0.75, /3 = 1. 




Figure 2. The left panel is the emef for a Pareto distribution with the parameters a = 7 
and A = 3 and the right panel is the one for an Exponential distribution with the parameter 
A = 2. 
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Figure 3. The erne/ for two Weibull distributions. The left panel concerns the one with 
the parameters /? = 1, r = 3.6 and the right panel concerns the one with the parameters 
P = 1.5, r = 0.2. 




Figure 4. The left panel is the emef for the Burr distribution with the parameters a = 0.5, 
A = 0.5, T = 5 and the right one is the emef for the Gomberz distribution with the parameters 
a = 1, A = 0.5. 






16 


MEAN EXCESS FUNCTION 



0 2000 4000 6000 



0.4 0.5 0.6 0.7 0.8 0.9 


Figure 5. The left panel is the emef for the Gamma distribution with the parameters 
a = 2, /3 = 0.001 and the right one is the emef for the Beta distribution with the parameters 
X = 7,/3 = 2. 




Figure 6. The left panel is the emef for the Lognormal distribution with mean /i = 0 and 
with variance cr^ = 1 and the right one is the emef for the Normal distribution with the 
same parameters. 






(n)u 8 (n)u i 
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Figure 7. Emef for Laplace distribution with the parameters p- = 0 ,(t = 1, r = 0.5 




u u 

Figure 8. Emef for a t Student with v = h degrees of freedom and skewing parameter 
p, = 0 (left panel) and for a Cauchy distribution with location parameter /r = 0 and scale 
parameter (5 = 1 (right panel). 
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5.2.2. Generalized hyperbolic models. Next, we consider the emefs for the Gh models. We need some defini¬ 
tions. The Lebesgue density function of the one dimensional Gh is given by 

(5.2) = a(A,a./3,5,M) ^ ^ X K^_i{ay/6^ + {x - fz)^) 

where 

(a^ — /3^)^ 

*■ ’ ^/^a^^~^l6^K\{Sy/a‘^ — /J^) 

is a norming constant to make the curve area equal to 1 and 

1 /I \ 

K\{x) = -J y^~^exp{^--x{y + y~^)jdy, (a; > 0) 

is the modified Bessel function of the third kind with index A. 


The dependence of the parameters X,a,P,S, and y is as follows: a > 0 determines the shape, 0 < |/3| < a 
the skewness, ya € M is a location parameter and <5 > 0 serves for scaling. The parameter A € M specifies the 
order K\ function Bessel that appears in the Gh density function and is used to obtain different subclasses 
of Gh distribution. 


In the following, we summarise the differents domains of possibilities for the parameters 


If 

if 

if 


A < 0, 

(5 > 0, 

|/3| < a 

A = 0, 

(5>0, 

l/^l < a 

A > 0, 

(5>0, 

l/^l < a 


An important Gh family aspect is that it embraces many special cases such that Hyperbolic (A = 1), 
Student-f(A < 0), Variance Gamma (A > 0), and the Normal Inverse Gaussian (NIG) (A = —0.5) distribu¬ 
tions. 

It nests also Generalized Inverse Gamma (GIG) distribution defined only by the three parameters A, a, and 
(3. An Inverse Gaussian (IG) distribution is a GIG distribution with A = —0.5 and a Gamma (T) distribution 
is also a GIG distribution with /3 = 0. 


It contains some limiting distributions such as Gauchy distribution with parameters y and S (obtained for 
A = —0.5 and a = /3 = 0). 


The Gaussian distribution with mean y and variance cr^ are obtained for A 


—0.5, for a, d —)■ oo and - —>■ cr^. 

’ ’ a 


The Skew-Student t with u degrees of freedom is obtained if a = |/3|, then u = — 2A > 0. 


The Student t distribution is obtained for a = (3 = 0, y = 0 and S = y/v. In the special case of hyperbolic 
distributions (A = 1), we obtain the skewed Laplace distribution for 5 = 0. 


All of these have been used to model financial returns and log-returns. 
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In TABLE we consider some specific Gh distributions with the superscript spe and limiting distributions 
with the superscript Im. The used parameters are specified and the emef figures corresponding to each choice 
are displayed. 


Distributions 


Parameters 

P 


Hyperbolique®^® 

1 

1.5 

-0.5 

0.75 

0.2 

t-Stud.®P*= 

-2 

10 "® 

0 

2 

0 

Asym. t-Stud.*P*= 

-1.278 

0.01186 

0.01186 

0.0766 

1.005 

-1.247 

0.0148 

-0.0147 

0.076 

1.005 

NIG**P*= 

-0.5 

8.03 

-1.37 

0.051 

0.0105 

-0.5 

7.6 

-1.24 

0.052 

0.0103 

Variance Gamma**^® 

2 

0.3 

0.1 

2 

0 

GiG‘>pe 

5 

3 

1 

- 

- 


IG Inverse Gaussian*"* 
Gamma(a, 
ir Inverse Gamma*"* 
Skew Laplace*”* 
Gaussian*"*(3, 0.3) 
Gauchy*"* (7,1) 


0.5 

0.5 

-0.5 

I 

-0.5 

-0.5 


1 

4.51012 
1.9 X 10"® 
1.1 
10 ® 

0 


0 1 

10 "® 

3.1 X 10"® 

0.1 0.001 

2 3 X 10® 

0 1 


2 

3 

7 


Figures 
Figure 
Figure 
Figure 


11 


Figure 12 


Figure 

Figure 

Figure 


Table 2. Specific and limiting GH distributions. 


5.2.3. Graphical test. We are now in a position to use the emef graphs already drawn as tools of goodness of 
fit. 


Emef for Normal Inverse Gaussian (NIG) and t-student-distributions are not monotonic function. They 
decrease and increase like for emef returns data. For this reason, we fit them to both monthly returns and 

FIGURE 


log-returns from Dow Jones data base (see FIGURE 17 
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FIGURE 


21 and FIGURE 23). 


Dow Jones data base consists of several compagnies like AXP(American Express compagny), GSGO(Gisco 
Systems), DAX, CAT, IBM and so one. Each one having 5 values : from opening {op) values to closing (cl) 
values , also minimum {min), maximum {max), and volume {vol) values. 


We select AXP and CSCO compagnies and we consider returns and log-returns for their values as showed in 
the TABLE Then we construct their emef plot and their fitted counterpart. 


Estimates parameters and the emef are given in table 
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Figure 9. Left panel : Emef for an hyperbolic distribution with the parameters A = 1, 
a = 1.5, /? = —0.5, 5 = 0.75, ^ = 0.2. Right panel : Emef for a t-student with the 
parameters A = —2, a = 10“®, /3 = 0, (5 = 2, /r = 0. 




Figure 10. The emef for two t-Student distribution. The left panel concerns the one with 
the parameters A = —1.278, a = 0.01186, [3 = 0.01186, 5 = 0.0766, ^ = 1.005 and the 
right panel concerns the one with the parameters A = —1.247, a = 0.0148, (3 = —0.0148, 
(5 = 0.07683, = 1.005. 
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Figure 11. The emef for two Normal Inverse Gaussian distributions. The left panel con¬ 
cerns the one with the parameters A = —0.5, a = 8.03, (3 = —1.37, 5 = 0.051, fi = 0.0105 
and the right panel concerns the one with the parameters A = —0.5, a = 7.6, /3 = —1.24, 
S = 0.052, n = 0.0103. 




Figure 12. Left panel : the emef for a Variance Gamma distribution with the parameters 
A = 2, a = 0.3, P = 0.1, (5 = 2, /i = 0. Right panel : the emef for a Generalized Inverse 
Gaussian distribution with the parameters A = 5, a = 3, ,d = l. 
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Figure 13. Left panel : the emef for Inverse Gaussian distribution with the parameters 
A = 0.5, a = 1, /? = 0, 5 = 1, /i = 0. Right panel : the Emef for a Gamma distribution with 
the parameters A = 0.5, a = 4.5 x 10^^, j3 = 10“^. 




Figure 14. Left panel : the emef for the Inverse Gamma distribution with the parameters 
A = —0.5, a = 1.9 X 10“®, /? = 3.1 x 10“^. Right panel : the emef for the Skew Laplace 
distribution with the parameters A = 1, a = 1.1, /3 = 0.1, 5 = 10“^, fi = 2 










(nlu i 
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u u 

Figure 15. The left panel concerns the emef for the Gaussian distribution with mean /i = 2 
and variance = 0.3. The right panel concerns the emef for the Cauchy distribution with 
location parameter n = 7 and scale parameter 5 = 1. 
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Compagnies 

Nature 

Values 

Emef plots 

Real emef 

Eitted emef 

AXP 

Returns 

op 

Figure 

16 

Figure 

17 


min 

Log-returns. 

max 

Figure 

18 

Figure 

19 


cl 

CSCO 

Returns 

min 

Figure 

20 

Figure 

21 


vol 

Log-returns 

op 

Figure 

22 

Figure 

23 


max 


Table 3. (Fitted) Emef for DAX and CSCO compagnies data. 


Comp 

Nature 

Values 

Ghyp estimates parameters 

Fit.Dist 

Figures 

A 

a 

/3 

A 

A 

AXP 

Returns 

op. 

-1.278 

0.01186 

0.0118 

0.0766 

1.005 

t-stud 

Fig. 

17 


min. 

-1.247 

0.0148 

-0.0148 

0.0768 

1.005 

t-stud 

Log-ret. 

max. 

-0.5 

8.03 

-1.37 

0.051 

0.0105 

NIG 

Fig. 

19 


cl 

-0.5 

7.6 

-1.24 

0.052 

0.0103 

NIG 

CSCO 

Returns 

min. 

-1.24 

0.0148 

-0.0148 

0.0768 

1 

t-Stud. 

Fig. 

21 


vol. 

-3.82 

4.22 

4.22 

0.613 

0.753 

t-Stud. 

Log-ret. 

op. 

-1.26 

0.83 

- 0.83 

0.07 

0 

t-Stud. 

Fig. 

23 


max. 

-1.32 

0.85 

-0.85 

0.076 

0 

t-Stud. 


Table 4. Emef for fitted Gh distributions to DAX and CSCO compagnies data. 


5.2.4. Commentaries. In view of figure 10 and figure 17 we can say that t studient distribution fits well 
opening and minimum values return for the American Express compagny AXP, whereas NIG distribution 
fits well maximum and closing log-returns values for the Cysco System compagny CSCO in view of FIGURE 
HD and FIGURE [m 
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u u 

Figure 16. Emef for AXP compagny (monthly data returns). The left panel concerns the 
opening values and the right one concerns the minimum values. 




Figure 17. A t-Student distribution is fitted to monthly data returns for AXP compagny 
(see figure [l6| . The left panel concerns a t distribution with the parameters A = —1.278, 
a = 0.01186, /? = 0.01186, S = 0.0766, and /i = 1.005 fitted to opening values. 

The right one concerns a t distribution with the parameters A = —1.247, a = 0.0148, 

P = —0.0148, 6 = 0.07683, /r = 1.005 fitted to minimum values. 
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Figure 18. Emef for AXP compagny (monthly data log-returns). The left panel concerns 
the maximum values and the right panel concerns the closing values. 




Figure 19. A NIG distribution is fitted to monthly data log-returns for AXP compagny 
(see Figure [T^. The left panel concerns the one with the parameters A = —0.5, a = 8.03, 

/3 = —1.37, i5 = 0.051,/i = 0.0105 fitted to maximum values. The right panel concerns the 
one with the parameters A = —0.5, a = 7.6, /? = —1.24, 5 = 0.052, /x = 0.0103 htted to closing 
values. 













MEAN EXCESS FUNCTION 


27 




Figure 20. Emef for CSCO compagny (data returns). The left panel concerns monthly 
minimum values and the right panel concerns monthly volum values. 




Figure 21. A t-student distribution is fitted to monthly data returns for CSCO compagny 
(see Figure [20|). The left panel concerns the one with the parameters A = —1.24, a = 0.014, 
/3 = —0.014, 6 = 0.076, ^ = 1 fitted to minimum values. The right panel concerns the one 
with the parameters A = —3.82, a = 4.22,/3 = 4.22, S = 0.613, /r = 0.753 fitted to volume 
values. 
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Figure 22. Emef for CSCO compagny (data log-returns). The left panel concerns monthly 
opening values and the right panel concerns monthly maximum values. 




Figure 23. A t-student distribution is fitted to monthly data log-returns for CSCO com¬ 
pagny (see Figure [2^. The left panel concerns the one with the parameters A = —1.26, 
a = 0.83, /3 = —0.83, 5 = 0.07, /i = 0 fitted to opening values. The right panel concerns 
the one with the parameters A = —1.32, a = 0.85,/3 = —0.85, S = 0.076, = 0 fitted to 

maximum values. 
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6. Conclusion 

In this paper we have established an asymptotic confidence bands for the mean excess function by using 
functional process approach. Then we applied these bands for fitting Gh distributions to Dowjones financial 
data. It is a known fact that these ones fit well financial data since they embrace major part of classic 
distributions. 

We remarked that Student and NIG distributions are good candidates for fitting returns and log-returns 
data showing their semi-heavy tails. 


7. Appendix 

7.1. Moment computations. 

Let Zi,Zn, n i.i.d centered random variables defined on the same probability space with common variance 
= Ki and common fourth moment E(Zj^) = K 2 > 0. We have 


E[T{n,u,S)]‘^ = e(^Zi + Z 2 + ... Zr^y 


= E 


k^l 




^E(Z|)+6 ^ E{ZlZ\ 

l<2<j<n 




(7.1) 

since, for distinct i, j, k and /, 

E{Z,Zf) = EiZ,Z]Zk) = EiZ,Z,ZkZi) = 0 , 

by using independence plus the fact that E{Zi) = 0. 

Using independence again, 

EiZfZ^) = E{Zf)E{Z^) = kI for z ^ j. 


We obtain ^ E{ZfZj) = since the number of possible couples {i,j) of integers such that 


l<2<j<n 


^ < i < j < n, is 


n\ n{n — 1) 


. Hence from ( |7.1[ ), we deduce that 
E[r(n, u, 5)] = nK 2 + 3n(n — 


□ 


7.2. Proofs of the uniform asymptotic consistency bounds. 

7.2.1. Talagrand bounds. We begin to recall the Talagrand bounds and a device of Einmahl and Mason on 
how to apply it. 

Before going any further, we recall that a class of measurable real valued functions J- is said to be a pointwise 
measurable class if there exists a countable subclass J-q of J- such as, for any function / in J-, we can find 
a sequence of functions {fm}m>o in for which fm{x) —>■ f(x) as m —>■ oo , a; S K. (See Example 2.3.4 in 

m- 

Further, let ^i, ^ 2 ) • • ■ be a sequence of independent Rademacher random variables independent of Xi,X 2 , ■ ■ ■, 
and Gm be the functional empirical process indexed by the class of functions X. 


The following inequality is essentially due to Talagrand (1994) (see [S]). 
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Inequality. Let be a pointwise measurable class of functions satisfying for some 
0<M<oo, ||/||oo<M, /GJ-, . 

Then for alH > 0 we have, 


max WVmGjnWj^ > + t) | 

l<m<n \ II T / J 

2 = 1 


(7.2) 


< 2( exp(—yl 2 + exp (—^2 t/M)), 


where CTjr = sup^gjr yar(/(X)) and Ai,A 2 are universal constants. 


And the lemma below of Einmahl and Mason is very helpful for obtaining bounds on this quantity, when 
the class has a polynomial covering number. 

Assume that there exists a finite valued measurable function G, called an envelope function, which satisfies 
for all cc G M, G{x) > sup |/(a;)|. We define for 0 < e < 1 


Nis,T) :=supiv(eVQ(G2),J-,dQ) 

Q ^ ^ 

where the supremum is taken over all probability measures Q on M for which 0 < (5(G^) := J {y)Q{dy) < oo 
and dg is the L 2 (Q)—metric. As usual, N{e,F,dQ) is the minimal number of balls {g : dQ{g,f) < e} of 
dg—radius e needed to cover T. Here is the device of Einmahl and Mason [9]. 


Lemma 3. (Einmahl - Mason [9]^ Let J- be a pointwise measurable class 
some constants /3 > 0, u > 0, G > 1, ct < 1/(8G) and function G as above, 


(A.1) 

(A.2) 

(A.3) 


E 


G^{X) 


< 


^2. 


N{e,T)<Ge-'', 0 < e < 1; 

o-Q := supE[/^(X)] < 

/e.T L J 


(A.4) sup ll/lloo < ^ A/ncr71og(/3V 1/g). 

/ejr 2Vi^ + l 

Then we have for some absolute constant A, 


of bounded functions such that for 
the following four conditions hold: 


(7.3) 


e|^6/(^.) 

2=1 


< A a/ uncr^ log(/3 V l/tr). 


7.2.2. APPLICATION. Put £u{x) = £{x)I(^x>u)^ with £(x) = 1 or £(x) = x, and let P = {£u,u G /}. 

P is pointwise measurable since it suffices to take Pq = {£u,u G / n Q}, where Q is the set of irrationnal 
numbers. 


Next G = max(|f'(uo)|, |^(ui)|) = M > 0 is an envelope of P since we have 

sup \£u{x)\ < \£{x)\ < max(|^(uo)|, I^(mi)|), y Uq < x < Ui- 

uei 

Remark that if £{x) = 1 then M = 1 and if £(x) = x then M = max(|ito|, li^il)- 
We have cr|r = sup Var{f{X)) < M^. 

So we may use Talagrand’s inequality : it remains to check points of Lemma 
Points (A.l) and (A.3) are obvious with (3 = M = a. 

To check (A.2), consider any probability Q on K. We get for (u,v) G P,u< v, 

d%{£u,£v) = Ji£u - £v f{x)dQ{x) < M'^Q{[u,v]). 


( 7 . 4 ) 
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By a classical result in probability in M, for any given 0 < £ < 1, we may cover [ug, mi] by at most 
m = sub-intervals such that s^]) < £^, i = 1,... ,to. 

([ai] stands for the smallest positive integer greater than or equal to x). 

Let C = {m + l)e, we have m < Ce~^. 


For any u G [ug, ui], there exists i G {1,..., m} such as Si_i < u < Si with Q{[u, s^]) < e^, so the corresponding 
iu G J- is such that 


<eM = £VQ(G2) 


from (7.4). 


To hnish m = N(^y^Q{G^, T, dq^ < Ce ^ and 


N{e,F) = supiV 
Q 


(e^/Q{^),F,dQ) 


< Ce-^. 


Now we take = max(2,max(|£(Mg)|, |t'(ui)|) = Mi. 

Finally for 

8M2 log Ml 


n > 


M2 


we have 

n 

(7.5) <G^v^, 

i=i ^ 

where Cjr = A Mi-^log Mi, since all the points of the Lemma are checked. 


Now we are going to apply the inequality (7.2) first for the class of functions 


= Uu(x) = gu{x),u G /}. 

In this case Mi = 2 since £{x) = 1, for any ug < a: < ui, and 


E = D„,i < where Cjr^ = 2A^\og2. 

t — 1 


Let £ > 0, ni > 2 log 2 and to such that 

ni 

(Remind that CTjt = 1.) 

Then 


exp ^and exp ^ — A 2 to^ < - and tg < y/^. 


i=l 


P| max WVmGmWj^i > Ai(E\\^^iguiXi) -h tg) | < e/2. 

So for u > ni, we arrive at 

IP(|IPn(5u) -Pa(5«)I 

As t^/y/n < yfn\lyfn < 1, we obtain 


< 


Al(7?n,l + to) 


,UG > 1 - e/2. 


thus 

(7.6) 


Ai + to) ^ A\Cjr^yfn -l- A\yfn ^ A\Cjr^ A\ ^ D\ 


|Pn(ff«) - Px(5«)l < ^ > 1 - £^/2 
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where Di = 2AAi^/\og2 + Ai. 


Let us use the same method, for the class of functions 

T 2 = {iuix) = fuix),u G /}. 

In this case Mi = max(2,max(|Mo|) |mi|)) since £{x) = x, for any uq < x < ui, and 

n 

E\\'^S_Ju{Xi) = Dn ,2 < CJr^^/n, where = AMuJlogMi. 

" i=i 

^ SM^logMi , , , 

Let 712 > -- and to such that 


exp { -< - and exp f — A 2 tol < ^ and tg < \A^- 

V 772 70 / 8 


Then 


n 

p| max \\^/m<Grn\\J ^2 > ^1 (eII II +^o)|<£/2- 

k l<ra<n \ N ^ 


So for 77 > 772, we deduce that 

P( \En{fu) - Ex{fu)\ < ^ uel) >l-e/2. 

As to/^/n < 'sjnil\/n < 1, we obtain 


thus 

(7.7) 


+ ^o) ^ + All /77 ^ AiCjr^ + Ai ^ Di 

~ 77 sjn sjn 


p(|P„(/„)-Px(/.)| < ^,77G/) >1-£/2, 


where = A Ai Mi/log Afi + Ai. 


Now we use again (2.3) : 

|e„(77) - e(u)| < |P„(/„) - Px(/'u)| X |P„( 5 „)r^ 
For UQ<u<ui^we get 


|Px(/«)| X 


|Pn(gn) - Px(g«)| 

|Pn(5^i)Px(5«)| 


0 < F(ui) 


\fn 


< Ex{gu) 


£1 

^/n 


< En{gu) < Exigu) + 


£l 


with a probability greater than (w.p.g.t) 1 — e/2 and thus 



so we obtain 

p(|Pn(/«)-Px(/«)| X |P„(5.)r' 

>P(|Pn(/.)-Px(/u)| <^,UGl) > l-e/ 2 , 

thus 

(7.8) P(^|Pn(/„)-Px(/„)| X |P„(5„)|-1 < ^ X (f(77i)-^)”',77G/) > 1 - e/2. 
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From the following inequalities : 

'‘n{9u)\~^ < (^F{ui) - , w.p.g.t l-c/2, 

\Vxi9u)\-^ <F{u,)-\ 

,\Fx{fu)\<nX\, 
and by the same manner, we obtain 


\Fn{guWx[9u)\ y/n 


xE|X| X (i^K)(FK)-^)) \u&I^ 


> 


]p( \Fn{9u) - Fx{9u)\ < ^ > 1 - e/2, 


thus 


(7.9) P |Px(/«)|x ' i^'V' ^ 7 - 7 ;' < ^ 

\Fn(9u)Fx{9u)\ y/n 


{gu) Px(5«)l ^ X E|X| X (^F{ui)iF{ui) - ^)) ') > 1 - e/2 


By combining (7.8) and (7.9), we obtain 


p(|Pn(5jr' X |P„(/„)-Px(/„)| > ^ X \ue?j 

ri9u)-Fxigu)\ ^ Di (F(wi)(F(mi)- 


p |Px(/.)|x ~7^T/' > ^ 

\Fn(gu)Fx[gu)\ y/n 


£ £ 

<- + -<£. 
~ 2 2 ~ 


This gives 


Er,,. . , „ 1 /„ DixE\X\ 


Kiu) - e{u)\ > ^) < £ where 

y/n F{ui) - Di/y/n 


D 2 + 


F{ui) 


Finally we conclude by : 


/ Eji E^i \ 

(e„(u)- X < e(y) < e„(u) + —F,u G I] > 1- £ 

y \/n Jn J 


for any £ > 0, any n > no with 
En = =- 


Fiui)-D 


i/y/n^ 


F(ui) 


Di = 2AAi y/\og 2 + Ai 
j and < D 2 = AAiMiy/\ogMi + Ai, 
Ml = max(2,max(|Mo|, |mi|)) 


□ 
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